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Abstract Confocal laser scanning microscopy was applied 
to visualize inside fibers stained with a fluorescent dye, and 
the fiber orientation distribution in the thickness direction 
was determined by image analysis from optically sliced 
images. Fiber orientation angle and anisotropy were deter¬ 
mined by the fast Fourier transform method, which holds 
advantages over the conventionally applied Hough trans¬ 
form with regard to evaluation of fiber width-weighted 
contribution and intrafiber segmental contribution. An 
orthogonally layered Japanese paper handmade by the 
flow-sheet forming method resulted in a clear change of 
fiber orientation angle from 0° to 90°. Machine-made wood- 
containing printing paper showed the highest anisotropy, 
while copier paper showed a low anisotropy. Thick papers 
like the wood-free paperboard and recycled packaging 
board showed that the maximum measured thickness was 
about 150 jum. The depth limit to detect fluorescence was 
considered to depend on the apparent density of the paper 
and the light absorption character of the fibers. 

Key words Confocal laser scanning microscopy • Fiber ori¬ 
entation • Handmade paper • Washi • Fourier transform 


Introduction 

In the papermaking industry, the use of innovative tech¬ 
nologies such as the lamination of paper with other novel 
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materials and embedding of electronic microchips will 
become more common in the near future. These new types 
of paper products will go outside the reach of the conven¬ 
tional planer sheet structure and must be dealt with as a 
three-dimensional bulk structure. Therefore, analytical 
techniques that consider the three-dimensional structure of 
paper are predicted to become increasingly important. For 
this purpose, several techniques have been established. 
Mechanical slicing of paper with a microtome to make cross 
sections is a reliable but time-consuming method. Thin 
cross-sectional slices made so can be observed with an 
optical microscope. Developmentally, consecutive smooth 
cross sections of epoxy-embedded coating layers after suc¬ 
cessive grinding and polishing were obtained and recon¬ 
structed as a three-dimensional representation to visualize 
pigment particles and pores. 1 Recently, the focused ion 
beam has emerged in place of the microtome to make intact 
cross sections. 2 

Fiber orientation is one of the major parameters that 
dominate the physical properties of paper. For example, 
curl occurs when paper is moistened. The curl axis is usually 
consistent with the fiber orientation axis. However, fiber 
orientation, per se, does not cause curl, but the two-sided 
nature or transverse variation of components or structural 
parameters including fiber orientation causes it. Thus, the 
way in which fiber orientation distributes in the transverse 
direction (z-direction) is worth measuring as one of the 
paper structural parameters to predict three-dimensional 
deformation. Actually, a flow of a fiber suspension induces 
fiber orientation on a forming screen, and its degree has a 
distribution. 3,4 

According to a survey of traditional handmade paper 
like washi (Japanese paper) and hanji (Korean paper), it 
was found that bast fibers tend to be oriented in the case of 
washi more remarkably than conventional machine-made 
papers. This is because a fiber flow on a slanted forming 
screen consisting of bamboo splints woven together with 
silk threads, corresponding to a “forming wire” in the 
modern papermaking industry, during dehydration results 
in evident fiber orientation as observed for the flow sheet 
forming. If this principle is used the other way, fiber orienta- 
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tion can be a physical basis to evaluate the effects of the 
motions of a mould (a wooden frame holding the forming 
screen) in the process of hand-making of cultural paper. 5 
Other than evaluation of hand-making processes, informa¬ 
tion on fiber orientation distribution of a repair paper in the 
thickness direction should be useful to predict deformation 
of the repair paper over time after repairing. 

Regardless of the sheet-forming tool (machine or hand), 
fiber orientation is one of the most important factors of 
paper properties and its distribution in the thickness direc¬ 
tion plays an important roll for deformation of the paper 
by both quick and slow environmental changes. To analyze 
the fiber orientation distribution, a representative conven¬ 
tional technique is the tape-stripping method. In this 
method, an adhesive tape is stuck on a paper surface; when 
the tape is ripped, a fiber layer with a certain thickness is 
removed together with the tape. The thickness and fiber 
orientation of the layer is measured to determine fiber 
orientation parameters as a function of the depth from 
the paper surface. Abe and Sakamoto 6 applied the tape¬ 
stripping method to separate thin layers from a whole sheet 
and measured Fraunhofer diffraction using a laser beam for 
fiber orientation. They confirmed that the angle-dependent 
distribution of optical density of the diffraction pattern 
recorded on the photographic film was in good agreement 
with fiber orientation distribution measured by another 
method and showed the characteristic fiber orientation dis¬ 
tributions in the thickness direction for a shaking wire 
machine and a high speed machine, respectively. As well, 
Erikkila et al. 7 applied the tape-stripping method to split a 
sheet to 10 to 12 individual layers to determine fiber orien¬ 
tation angle and the anisotropy for each layer. The method 
for establishing the magnitude and direction of fibers in 
images taken with a video camera was image analysis using 
a Sobel operator-based procedure. Although such local 
filtering operator-based procedures are excellent to ac¬ 
complish a quick operation, thick fibers are not weighted 
depending on their fiber width because only edges are 
detected sensitively by those filtering operators. 

In the tape-stripping method, fiber breakage and pullout 
is inevitable and there is no guarantee that the layer thick¬ 
ness is uniform. A novel technique employed as an alter¬ 
native in the present work is confocal laser scanning 
microscopy. The confocal system provides optical (non¬ 
mechanical) slicing, which solves the problems typical 
of mechanical layering, although sufficient precision is 
restricted within reach of the emitted light from fibers. 
Sliced layer images are obtained directly and are easily 


subjected to image analysis. As an image analysis tech¬ 
nique, Xu et al. 8 applied the Hough transform to detect fiber 
edges in pioneering work to acquire the index and angle of 
fiber orientation from confocal images. However, this 
popular method includes several problems in determination 
of fiber orientation. Only two edge lines are detected by the 
Hough transform for one fiber whether the fiber is thick or 
thin as well as the Sobel operator. In addition, curved fibers 
tend to be excluded from orientation calculation when the 
segments constituting one fiber are short, although each 
segment contributes to its tangential direction depending 
on its length. Therefore, in this work, a fast Fourier trans¬ 
form (FFT)-based algorithm was applied to determine fiber 
orientation parameters. The validity of the FFT method to 
evaluate fiber orientation was proved by Yuhara et al. 9 with 
the result obtained by applying it to a soft X-ray micro¬ 
diagram of a paper sheet. This technique was extended to 
the three-dimensional fiber orientation of paper in our 
research. 


Experimental 

Samples 

Five paper samples were used. Table 1 lists these samples 
with basis weight, thickness, and apparent density. The 
orthogonally layered Japanese paper was manufactured by 
Mr. Makoto Numata, a craftsman of Mino City, Gifu Pre¬ 
fecture. He produces the paper called sai by layering two 
thin sheet plies orthogonally to give one finished sheet of 
paper with improved dimensional stability. The other paper 
samples are commercially available machine-made papers. 

Two test pieces measuring about 20 x 20 mm were taken 
from each sample. They were soaked in a 0.02% ethanolic 
solution of acridine orange, which is a fluorescent dye to 
stain pulp fibers. Each test piece was placed on a slide glass 
and covered with a cover glass. Then a drop of aniline was 
put on its corner or edge to allow it to penetrate the sample 
spontaneously to the other end, leaving as few air bubbles 
as possible. As an impregnation liquid in the case of pulp 
fibers, aniline was selected because its refractive index is 
1.58, which is close to that of pulp fibers ranging between 
1.57 and 1.60. A light beam does not refract at an interface 
between substances with equal refractive indices, but passes 
straight. Even a porous material like paper appears trans¬ 
parent when such a liquid fills its inside pores. Note that 


Table 1 . Paper samples used in this study 


Paper 

Basis weight (g/m 2 ) 

Thickness (pm) 

Apparent density (kg/m 3 ) 

Orthogonally layered Japanese paper 

24 

55 

440 

Wood-containing D-grade printing paper 

51 

87 

590 

Copier paper, 100% recycled 

64 

96 

660 

Wood-free paperboard 

155 

167 

930 

Packaging paperboard, 100% recycled 

365 

453 

810 
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Fig. 1 . Acquisition of three-dimensional image by confocal laser scan¬ 
ning microscopy (CLSM) 

aniline is toxic and care must be taken not to touch it or 
inhale its vapor. 

Microscope conditions 

A confocal laser scanning microscope (CLSM) LSM 510 
with an upright body (Axioplan 2, Carl Zeiss, Germany) 
was used. As Fig. 1 shows, functions characteristic of this 
type of optical microscope to obtain three-dimensional 
(3D) images are optical slicing in the transverse direction 
by the confocal system and the intense power of laser beams 
to compensate for the resultant insufficient illumination. 
The term “confocal” means collecting light exclusively from 
a single plane in focus with a pinhole that eliminates light 
reflected from planes other than the focal plane. Laser 
beams fall on the front of a sample, and then reflect and 
enter the detector to provide a fluorescence image. The 
lenses used were a lOx objective and a lOx ocular to observe 
the sample at lOOx magnification. 

Sample observation used laser irradiation at 488 nm and 
optical filtering to restrict detection of light to fluorescence 
at about 512 nm that was emitted by fibers stained with 
acridine orange. The dimensions of the acquired images 
were 0.92 x 0.92 mm at a resolution of 512 x 512 pixels. 
Confocal plane thickness was 5 jum for paper and 8 jum for 
paperboard. The number of optical slices was 21 to 29 for 
one location and three to five locations were measured to 
average fiber orientation parameters at each depth. 

Paper thickness measured by this CLSM method is com¬ 
parable with that measured by the micrometer method, 
because the clamping pressure that develops between a 
cover glass and a slide glass is of the same order as that 
exerted by the micrometer platen. The clamping pressure 
P for the CLSM method depends on the surface tension y 
of the liquid and the distance d between the cover and slide 
glasses. P is defined as 2 y/d. y of aniline is 42.5 mN/m. If d 
is assumed to be 100 jum , P is calculated to be 85 kPa, which 
is equivalent to P specified in the standard 10 for measuring 
paper thickness (50 or 100 kPa). 


Image analysis to determine fiber orientation properties 

Acquired optically sliced layer images were analyzed to 
determine the angle and anisotropy with regard to fiber 
orientation. Two-dimensional FFT was applied to deter¬ 
mine those parameters. Fiber orientation can, in principle, 
be defined as an angular distribution of a power spec¬ 
trum as a result of the integration q(0) of a power spectrum 
P(r, 6) in polar coordinates, as shown in Eq. 1, where r 
is the radius and 6 is the polar angle. In the actual calcula¬ 
tion procedure adopted in this work, the images were 
binarized using a dynamic threshold method of a simple 
moving average. This was performed to correct the shading 
of images taken under nonuniform illumination that would 
frequently obscure the measurement of true fiber 
orientation. 

N/2 

q(0) = ^P(r,e) (1) 

r =0 

A two-dimensional Fourier transform inevitably pro¬ 
duces irregular large powers only in the x and y axes in the 
frequency domain when the gray levels of the left and right 
edges and the top and bottom edges of the image differ 
greatly, because FFT assumes infinite repetitions of the 
image in those directions. This binarization process effi¬ 
ciently reduces this abnormally large power by assigning an 
almost identical number of pixels to black or white. In addi¬ 
tion, this method of binarization can extract the edges of 
out-of-focus fibers located below in-focus top fibers even in 
completely out-of-focus images. The FFT was computed 
with these binary images to obtain the power spectra. 

When mean amplitude values are determined for each 
value of angle 6 , each discrete frequency value of the 
Fourier coefficients in the X-Y coordinates cannot be easily 
converted to the exact corresponding position in the polar 
coordinates. If one draws the radius of a given angle from 
the pole in the frequency domain, it would pass through 
only a few discrete positions. The interpolation was applied 
using four adjacent amplitude values. The procedure is 
detailed in the literature and is summarized in the follow¬ 
ing way. In the polar coordinates, the amplitude of the 
Fourier coefficient was added in a radial direction from the 
pole and its mean A{6) was determined for each value of 
angle 6 ranging from 0° to less than 180°. The closed curve 
combining A(0) and A(0 + 180) = (AO) was approximated to 
an ellipse using FFT followed by inverse FFT after low-pass 
filtering. The maximum value a , value b in the direction 
crosswise to a , and the angle a , that is, 0 when A(0) = b , 
were adopted as the lengths of the major axis, minor axis, 
and the orientation angle for elliptic parameters, respec¬ 
tively. Consequently, the anisotropy and the fiber orienta¬ 
tion degree are defined as alb and a , respectively. 

For one particular sample, representative values of ori¬ 
entation angle and anisotropy should be determined. Sliced 
layer images at the same depth of the same sample were 
collected from several locations within the same sample 
sheet. The mean orientation angle and anisotropy are 
defined as the means of several orientation angle and anisot¬ 
ropy values, respectively, calculated from every image at 
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the same depth. However, there was a problem in calcu¬ 
lating the mean orientation angle. Orientation angle is 
expressed a value between 0° and less than 180°. The angles 
of 0° and slightly less than 180° represent the almost the 
same orientation angle, but the mean angle is calculated to 
be about 90°. In order to avoid this misleading result, the 
overall calculation was devised. The overall orientation 
angle and anisotropy is defined as a single value calculated 
overall from accumulated amplitude values of Fourier coef¬ 
ficients at the same angle of every image at the same depth. 
The difference in the orientation angle and anisotropy 
made by the two calculation methods is shown in the 
results. 


Depth from surface 


Original slice image 


Binary slice image 


(a) 25 pm 


(b) 40 pm 


(c) 55 pm 


(d) 70 pm 


(e) 85 pm 



Fig. 2a-e. Original and binary confocal slice images of fibers in cross- 
lapped Japanese paper sai 


Results and discussion 

Figure 2 shows a set of optically sliced layer images for 
every 15 jum in depth for the orthogonally layered Japanese 
paper from the top side of the first ply to the screen side of 
the second ply. In this context, the screen side corresponds 
to the wire side in the modern papermaking industry. The 
images shown are actually every third image from a series 
of 25 images taken for slices every 5 jum apart. The fibers in 
the top two images (Fig. 2a, b) are oriented in the vertical 
direction of the image while those in the bottom (screen 
side) two images (Fig. 2d, e) are oriented in the horizontal 
direction. The central image (Fig. 2c) appears to include 
equal numbers of fibers oriented in both the vertical and 
horizontal directions. The orthogonally layered planes 
appear to be located near the image shown in Fig. 2c. The 
sliced layer images above that shown in Fig. 2a or below 
that shown in Fig. 2e gave very faint fiber shapes, probably 
because leak fluorescence emitted by fibers in adjacent 
planes was detected, even though there was practically no 
fibers in those actual confocal slices. The thickness of this 
sample was measured to be 55 jum with a micrometer, which 
roughly corresponds to the difference in depth between 
images shown in Fig 2a and Fig. 2e. For all of the sliced 
layer images, fiber orientation angle and anisotropy were 
calculated using FFT-based image analysis. 

Figure 3 shows a calculated result of fiber orientation for 
one location of the orthogonally layered Japanese paper. 
The fiber orientation angle is consistently around 90° for 
the top half (55 jum thick), and around 0° for the bottom 
half (65 jum thick), which is consistent with the visual per¬ 
ception of the sliced layer images in Fig. 2. The thickness of 
this orthogonally layered paper is 55 jum. Thus, outside this 
thickness, light reflected from the top or bottom surfaces is 



Depth from top, pm 


Fig. 3. Orientation angle (diamonds, left axis) and anisotropy (circles, 
right axis) of fiber orientation for one location of cross-lapped Japanese 
paper. Dotted lines marked a to e correspond to depths of slice images 
shown in Fig. 2 
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Depth from top, jum 


Fig. 4. Mean orientation angle ( diamonds, left axis) and mean anisot¬ 
ropy ( circles , right axis) of fiber orientation for cross-lapped Japanese 
paper. Bars show 95% confidence intervals 



Depth from top, pm 


Fig. 5. Overall angle ( diamonds, left axis) and anisotropy (squares, 
right axis) of fiber orientation for cross-lapped Japanese paper 


considered to be detected faintly. The confocal mechanism 
does not guarantee complete limitation of light from out- 
of-focus planes. Then indistinct virtual images are detected 
about 20 jum above and below the paper surface layers. 
Therefore, the paper is measured to be about 40 jum thicker 
by the CLSM method than the thickness measured by the 
micrometer method. Rough papers like Japanese paper 
with fibers rising from the surface are measured to be about 
10 jum thicker. The anisotropy is sensitive to the sheet struc¬ 
ture change of orthogonal layering. At 55 jum from the top, 
the anisotropy was lower than 1.1. With this calculation 
process, anisotropies higher than 1.15 suggest strong fiber 
orientation while those lower than 1.10 imply almost random 


Depth from surface 


Original slice image 


Binary slice image 


(a) 25 pm 


(b) 40 pm 


(c) 55 pm 


(d) 70 pm 


(e) 85 pm 



Fig. 6a-e. Original and binary confocal slice images of fibers in wood- 
containing paper 

orientation. Thus, fiber orientation near the lapping plane 
was estimated to be random due to the mixture of fibers 
aligned in cross directions. This sheet has two peaks of 
anisotropy in the centers of the two respective plies. 

Figure 4 shows the mean fiber orientation angle and 
anisotropy calculated from ten sets of sliced images for the 
orthogonally layered Japanese paper. Each set of images 
was shifted in the depth direction so that the angle change 
plane was between 35 and 40 jum in depth. There was no 
difficulty during this arrangement because every set of the 
sliced layer images had the plane exhibiting a clear orienta¬ 
tion angle change. Figure 4 confirms that this orthogonally 
layered paper has the crosswise structure over the whole 
area and that the CLSM technique is feasible to determine 
the internal fiber orientation. 
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Depth from top, pm 


Fig. 7. Overall angle ( diamonds, left axis ) and anisotropy (squares, 
right axis) of fiber orientation for wood-containing paper 


Figure 5 shows an overall fiber orientation angle and 
anisotropy of the orthogonally layered Japanese paper. The 
data for Fig. 4 are mean values of angle and anisotropy of 
all images at each depth. Therefore, 95% confidence inter¬ 
vals can be calculated, as shown in Fig. 4, from all the 
values. However, the data for Fig. 5 are single values of 
angle and anisotropy, and it is impossible to calculate 95 % 
confidence intervals. This overall calculation means the 
large-area sampling that covers the sum of the imaged areas 
subjected to the analysis and permits determination of 
the true orientation angle. Overall anisotropies tend to be 
lower than the corresponding mean anisotropies. Overall 
anisotropies higher than 1.10 suggest strong fiber orien¬ 
tation while those lower than 1.05 imply almost random 
orientation. 

Figure 6 shows a set of sliced layer images for wood- 
containing printing paper. Among machine-made paper 
samples used in this work, this wood-containing paper is 
likely to have the highest anisotropy. The fibers are observed 
to be orientated clearly in the machine direction, namely, 
the horizontal direction in the image. Figure 7 shows 
the overall angle and anisotropy of fiber orientation for the 
sample. The highest anisotropy reaches 1.30 across the 
thickness. The sliced layer with this anisotropy peak occurs 
toward the top surface, which is probably the wire side, 
considering that the wire side of paper is formed earlier 
than the top side, immediately after fibers contact the wire, 
causing more remarkable fiber orientation. 

Figure 8 shows a set of sliced layer images for copier 
paper. Generally, copier paper is designed for less fiber 
orientation to prevent wrinkles that may occur around 
heating rolls that are installed in copying machines to fix 
toner. For the sliced layer images, fiber orientation is not 
discernible by eye. Figure 9 shows that the anisotropy is 


Depth from surface 


(a) 25 pm 


(b) 40 pm 


(c) 55 pm 


(d) 70 pm 


(e) 85 pm 


Original slice image 



Binary slice image 










Fig. 8a-e. Original and binary confocal sliced layer images of copier 
paper 


lower for the copier paper than any other paper samples 
used here. 

Figures 10 and 11 show overall angle and anisotropy of 
fiber orientation for the wood-free paperboard and recycled 
packaging board, respectively. For the wood-free paper- 
board, the possible range of measured thickness was 152 jum 
at the highest, although the actual thickness measured by 
micrometer was 167 jum. Reflected light photographed in 
the sliced layer images near the last image at 152 jum in 
depth is so faint that sliced layer images below this layer 
would not provide accurate analysis. Paper with very high 
density like this wood-free paperboard absorbs and scatters 
light reflected from a deep layer to a great extent. The 
results for the two paperboard samples suggest that internal 
structure analysis by CLSM is applicable for paper or paper- 
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Depth from top, pm 



Fig. 9. Overall angle ( diamonds , left axis) and anisotropy ( diamonds , 
right axis) of fiber orientation for copier paper 


Depth from top, pm 

Fig. 11. Overall angle ( diamonds , left axis) and anisotropy ( diamonds , 
right axis) of fiber orientation for paperboard for paper container 
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Depth from top, pm 


Fig. 10. Overall angle ( diamonds , left axis) and anisotropy ( diamonds , 
right axis) of fiber orientation for wood-free paperboard 


board to a maximum depth of 150 jum , although the actual 
depth depends on the apparent density. The packaging 
board had a higher maximum value of anisotropy than the 
wood-free paperboard. Below the sliced layer at 64 jum 
deep, the anisotropy decreased with the depth for the same 
reason; namely, an insufficient amount of reflected light. In 
particular, dark gray fibers used in recycled paperboards 
tend to absorb emitted fluorescent light. 


stained with acridine orange. The fiber orientation distribu¬ 
tion in the thickness direction was determined by applying 
Fourier transform image analysis to every sliced layer image 
of paper. Fiber orientation angle and anisotropy were 
determined by the FFT method, which holds advantages 
over the Hough transform with regard to evaluation of fiber 
width-weighted contribution and intrafiber segmental con¬ 
tribution. Japanese handmade paper is known to have a 
strong orientation because flow of the fiber suspension ori¬ 
entates fibers during dehydration in flow-sheet forming. An 
orthogonally layered Japanese paper, a model of a structure 
with a clear change of fiber orientation at the orthogonally 
layered plane, resulted in a clear change of fiber orientation 
angle from 0° to 90° and two peaks of anisotropy in the two 
plies, respectively. This result demonstrates the validity of 
the CLSM method. For machine-made papers, the wood- 
containing printing paper showed the highest anisotropy, 
while the copier paper showed a low anisotropy presumably 
because copier paper is designed to be free of wrinkles 
during copying. Thick papers like the wood-free paper- 
board and recycled packaging board showed the possible 
range of measured thickness was about 150 jum. Fluores¬ 
cence emitted from fibers deeper than 150 jum is too faint 
to be used for FFT-based image analysis, although the 
depth limit depends on the apparent density of the paper 
and the light absorption character of the fibers. 

Acknowledgments The authors thank Mr. Makoto Numata for pre¬ 
paring orthogonally layered Japanese paper sai. 


Conclusions 

Confocal laser scanning microscopy is a useful technique to 
examine the internal structure of paper or paperboard when 
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